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ABSTRACT 

In January 2009, the 2.1-sec anomalous X-ray pulsar IE 1547.0—5408 evoked intense 
burst activity. A follow-up Suzaku observati on on January 28 re corded enhanced per- 
sistent emission both in soft and hard X-rays (|Enoto et al.ll2010~b[) . Through re-analysis 
of the same Suzaku data, 18 short bursts were identified in the X-ray events recorded 
by the Hard X-ray Detector (HXD) and the X-ray Imaging Spectrometer (XIS). Their 
spectral peaks appear in the HXD-PIN band, and their 10-70 keV X-ray fluenccs range 
from ~2x 10 -9 erg cm -2 to 10 -7 erg cm -2 . Thus, the 18 events define a significantly 
weaker burst sample than was ever obtained, ^10 _8 -10 -4 erg cm -2 . In the ~0.8 to 
^300 keV band, the spectra of the three brightest bursts can be represented success- 
fully by a two-blackbody model, or a few alternative ones. A spectrum constructed 
by stacking 13 weaker short bursts with fluences in the range (0.2-2) xlO -8 erg s" 1 
is less curved, and its ratio to the persistent emission spectrum becomes constant at 
~170 above ^8 keV. As a result, the two-blackbody model was able to reproduce 
the stacked weaker-burst spectrum only after adding a power-law model, of which 
the photon index is fixed at 1.54 as measured is the persistent spectrum. These re- 
sults imply a possibility that the spectrum composition employing an optically-thick 
component and a hard power-law component can describe wide-band spectra of both 
the persistent and weak-burst emissions, despite a difference of their fluxes by two 
orders of magnitude. Based on the spectral similarity, a possible connection between 
the unresolved short bursts and the persistent emission is discussed. 

Key words: stars: pulsars: general - pulsars: individual (IE 1547.0—5408, 
SGR J1550-5418, PSR J1550-5418, G327.24-013) 



1 INTRODUCTION 

Magnetars are a peculiar sub c lass of isolated neutron star s 
jDuncan fc Thompson! 1 19921 : iThompson fc Duncan! ll995T ) 
with evidence for ultra-strong magnetic fields, mainly emit- 
ting in the X-ray frequency. Located along the Galactic 
plane and in the Magellanic Clouds, ~9 soft gamma re- 
peaters (SGRs) and ~12 anomalous X-ray pulsars (AXPs) 
are currentl y believed to be of the magnetar class (fo r 
reviews, see IWoods fc Thompson! 120061 : iMereghettil l200St) . 
X-ray pulsations have been observed from them in a pe- 
riod range of P = 2-12 s, indicating slowly rotating iso- 
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lated pulsars. Together with their large period derivatives 
P ~ 1CP 12 — 10~ 10 s s _1 , the magnetic-dipole-radiation 
approximation indicates these objects to have ultra-strong 

magnetic fields of B = 3.2 x 10 19 ^/p~P G > 3.3 x 10 13 G. 
Thus, it is widely believed that the peculiar characteristics 
of SGRs and AXPs originate from their extremely strong 
magnetic fields. 

The X-ray radiation from magnetars emerges on a large 
variety of time scales and intensities. One form is "persis- 
tent" X-ray emission with a typical luminosity of L x ~ 10 35 
erg s _1 , which is usually stable over long periods of time 
(~ a few month or longer). A typical persistent X-ray lu- 
minosity of magnetars exceeds by 1-2 orders of magnitude 
that available from their rotational energy losses, ~ ifj 33-34 
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erg s _1 . Without evidence for mass accretion, the emission 
is therefore considered to be powered by dissipation of the 
magnetic energies. The broad-band persistent X-ray spec- 
trum is generally composed of blackbody-like soft X-ray 
component with a temperature of ^0.3-0 .5 keV (|Mereghettil 
l200Sl : lRea et al.ll2007l ; lEnoto et alj|201ll ). and a power-law- 
lik e hard X-ray component with a photon i n dex of V ~ 
1 jKuiper et al.1 12009 ; Iden Hartog et all 120081 ; lEnoto et al.1 
l2010d ). While the soft component can be considered as 
optically-thick thermal radiation from the stellar surface, 
the hard component, spanning from ~10 keV to ~200 keV 
or higher, is considered to emerge through a different and yet 
unidentified process (e.g.. [Thompson fc Beloborodovl 2005 
3nl 



2007 



Beloborodov fc ThompsonT 20071; Baring fc Harding 
Fernandez fc Thompsonll2007l : iHevl fc Hernauisdl2005l ) 

Another form of X-ray radiation from magnetars is spo- 
radic emission of bursts with a typical duration from ~0.1 
seconds to a few hundred seconds. These bursts are phe- 
nomenologically classified into three types: quite rare "giant 
flares" with L x > 10 45 erg s -1 lasting about a few hundred 
secon ds (jMazets et aI.lll979l : lFeroci et al1l200ll ; lHurlev et al.1 
120051 ). "intermediate fl ares" with L x ~ 10 42 - 10 43 erg s" 1 
lasting a few seconds (Kouveliotou et al 1 l200ll : [Ol ive et al.l 
|2004 llsrael et alj|2008r ). and much" more frequently occur- 



-10 41 ergs 



with ~0.1- 



ring "short bur sts" with L x 

sec durations (jNakagawa et al ] l2007bl : llsrael et al.l l200&t ). 
Thus, these explosive events often show luminosities exceed- 
ing the Eddington limit for a neutron star of 1.4M© (M© 
being the solar mass), Z/Edd ~ 1.8 x 10 38 erg s _1 , presumably 
due to suppressi on of the electron scattering cross sections in 
the strong field (|Paczvnskilll992[ ). The mechanisms of these 
bursts are thought to be related to re arrangement of the 
magnetic fields invoking reconnection (|Lvutikov|[2003T) . or 
motion and fracturing of the neutron star crust, i.e., star- 
quake [(Thompson et al-llioolj ). 

The persistent X-ray intensity of a magnetar sometimes 
increases by 1-2 orders of magnitude with unpredictable 
timing. Such an transient enhancement lasts typically for 
a few months, including its gradual decay. They are of- 
ten accompanied, at its early phase, by a burst activity, 
which can even lead to discoveries of new magnetars. So far, 
such burst-active states, or outbursts ; were observed from 
some magnetars; XTE J1810-197 (IGotthelf et al.l I200I 



20051; 



Ibrahim et al] 12004 llsrael et al.l I200I IGotthelf fc Halpernl 

CXOU J164710.2-455216 



Bernardini et al.l 



Muno et al.l 



Enoto et al.l 
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Figure 1. Long-term X-ray monitoring of IE 1547.0—5408 with 
the Swift/XHT, produced from the public Swift data archive 
through a standard analyses procedure, fitted by an absorbed 
blackbody spectrum with the absorption left free, (a) The ab- 
sorbed 2—10 keV X-ray flux is shown since 2007 May. The Suzaku 
observation is shown in the star (green), (b) An expanded view 
of panel (a) during the 2009 outburst. The decay is shown as 
a double logarithmic plot sinc e the first short burst (2009-01-22 
01:32:41; iGronwall et aDl2009h . 



Rea et al.l 



20091; 



l2003l ; lNakagawa|[2007al ). Such a possibility has been ex- 
amined using a cumulative number-intensity distributio n 
of short bursts (jGotz et al.1 120061 : iNakagawa et al.l l2007bh . 
However, the observational information has so far remained 
insufficient to evaluate this possibility, since the studied 
short bursts are so bright (with fluence > 10~ 7 ergs cm -2 ) 
and infrequent that their time-averaged flux is much lower 
than that of the persistent emission. Therefore, it is inter- 
esting to examine, from observations of activated magne- 
tars, whether weaker short bursts become similar in spectral 
shape to the pers istent X-ray emission, as recently found in 
SGR 0501+4516 jNakagawa et al.ll201ll ). 

In the present paper, we studied the activated magne- 
tar IE 1547.0 —5408. Following its di scovery by the Einstein 
Observatory (lLamb fc Markert Il98lf) in 1980 and confirma- 



2007 : llsrael et al l 120071), SGR 05 1+451 6 tion by ASCA (|Sugizaki et al 



Enoto et al. 2010a) 



200 lh . this object was rec- 



andlE 1547.0-5408 (|Mereghetti et al.l |2009| ; lEnoto et al.l 



l2010bl ). More recently, there have been accumulating reports 
on such activities from sources with much weaker dipole 
fields (< 4.4x10" G) ; SGR 0418+5729 (Ivan der Horst et al l 
201fj:lRea et al.ll2010l h SGR 1833-0832 llGogiis et al Il2010l : 
Esposito et al.ll201lh . Swi ft J1822. 3-1606 jRea et al.ll2012T ). 
and Swift J1834.9-0846 (|Kargaltsev et alj|2012h . 

The enhanced persistent emission and the burst activity 
have been simultaneously observed in many activated mag- 
netars. However, it is not yet clear how these two emission 
forms are physically related with each other in the postu- 
lated dissipation process of the magnetic energy. One inter- 
esting possibility is that the persistent emission is composed 
of a large number o f small bursts that are not individu- 
ally detectable (e.g., iThompson fc Duncan! [l996l ; iLvutikovl 



ognized as a magnetar candidate, loca ted at the center of 
the s upernova remnant G327.24— 013 (|Gelfand fc Gaenslerl 
l2007h . based on its X-ray spectrum from XMM-Newton 
and Chandra observations. This was followed by a discov- 
ery of its radio pu lsations at ~2.07 s (PSR J1550-5418; 
ICamilo et al1l2007f), and o f X -ray pulsations at the same 
period 1 Halpern et alJlioOfet ) . Together with a measured pe - 
riod derivative, P ~ 2.3 x 10" 11 s s _1 (|Camilo et alj|2007h . 
its characteristic age and surface magnetic field intensity 
were estimated to be 1.4 kyr and 2.2 x 10 14 G, respectively; 
both parameters fall in the magnetar regime. At present, 
IE 1547.0—5408 is one of the fastest rotating objects among 
the known mangetars. 

Figure [T] shows a long-term monitoring of the 2- 
10 keV persistent emission from IE 1547.0—5408. In 
2007, a small enhancement was observed in its persis- 
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Figure 2. Background-inclusive count rates of IE 1547.0—5408 
during the present Suzaku observation performed on 2009 January 
28. From top to bottom, panels refer to those obtained with XIS1 
plus XIS3 (2-10 keV), XIS0 (2-10 kcV), HXD-PIN (10-70 keV), 
and HXD-GSO (50-150 kcV). Data are binned into 1 sec, except 
for 2 sec of panel (a). 



tent luminosity (|Halpern et alj 120081 ). In 2008 October, 
short bursts and decayi ng persistent sof t X-rays have 
been monitored by Swift |lsrael et al.ll20icil ). Three month 
later, a much stronger burst activity, to be dealt with 
in the present paper, was detected on 2009 January 22 
(UT). A series of intense short bursts were recorded 



satellites; Swift dGronwall et al. 


120091). IN- 


(|Savchenko et al.l 2009), Suzaku ( 


Terada et alj 



von Kienlin fc Connaughtonl 12009)1, RHESSI (iBellm et all 
20091 ), and Konus- Wind jGolenetskii et all 120091 ). Some 



bursts were very bright with fiuences above ~ 10 erg 
cm " 2 (|Mereghetti et al.ll2009l ), and a 150-s-long enhanced 
quasi-persistent emission was also recorded o n one occa- 
sion by the Fermi/GBM (|Kaneko et all I2010T ). Based on 
these SGR-like act ivities, IE 1547.0—5408 was also named 
SGR J1550-5418 (|Kouveliotou et al.ll2009l ). 

As shown in Figure [1] the early burst activity was ac- 
companied by a clear enhancement of the persistent emis- 
sion, which decayed on a timescale of a few month. In Enoto 
et al. (2010b; hereafter Paper I), we studied persistent emis- 
sion of this object using the Suzaku data acquired on 2009 
January 28-29, or ~7 days after the burst onset (Figure [T]), 
when the 2-10 keV persistent flux became by 1-2 orders 
of magnitude higher than that during less active states in 
2006 and 2007. On that occasion, a persistent hard X-ray 
component was detected for the first time from this source 
at least up to ~110 keV. The acquired broadband spec- 
trum in 0.7-114 keV was reproduced by an absorbed black- 
body (BB) emission with a temperature of 0.65 ± 0.02 keV, 
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Figure 3. (al) A At distribution of the 10-70 keV HXD- 
PIN events. Dashed-dottcd line indicates the best fit exponential 
model with the average event rate of A = 0.778±0.001 count s . 
(a2) An expanded view of panel (al) in the 0-0.3 sec range, (bl) 
Distribution of the consecutive burst-like counts calculated using 
panel (a). (b2) An expanded view of panel (bl). Dashed-dottcd 
line represents the predicted chance coincidence probability of 
ordinary events, (O.OSG)^ 138 ' . 



plus a prominent hard power-law (PL) with a photon index 
of r per ~ 1.5. The enhanced persistent emission was also 
studied by Chandra, XMM-Newton, Sw i ft, and INTEGRAL 
llBernardini et al.l l201ll: iNg et alJl201ll: Kuiper et al.l [20091 : 



Iden Hartog et a l. 2009; B ernardini et al.ll201ll ) 



In this paper, we present spectral studies of short bursts 
detected with Suzaku in the 2009 observation. The low back- 
ground of Suzaku, especially in the hard X-ray band, allowed 
us to detect much weaker short bursts than observed so far 
from any magnetars. Unless otherwise specified, we show all 
uncertainties at the 68% (la) confidence level, 



2 OBSERVATION AND DATA REDUCTION 

2.1 Observation and Data Screening 

We utilized the same Suzaku data (OBSID 903006010) of 
IE 1547.0—5408 as used in Paper I, although short burst 
events were eliminated therein. The data were acquired over 
~1 day from 2009 January 28 21:34 (UT) to January 21:32 
(UT). As described in Paper I, two (XIS1 an d XIS3) among 
the t hree X-ray Imaging Spectrometer (XIS: lKovama et al.l 
sensors were operated with a 1/4 window mode plus 
burst option, repeating a 0.5-sec exposure and a 1.5-s artifi- 
cial dead-time. The other sensor, XIS0, was operated in the 
timing mode (P-sum mode) to attain a high time resolution 
of ~7.8 ms together with one dimension al position infor- 
matio n. The Hard X-ray Detector (HXD: iTakahashi et al.l 
120071) . composed of silicon PIN diodes (HXD-PIN; 10-70 
keV) and GSO scintillators (HXD-GSO; ~50-600 keV), was 
operated in the normal mode with 61 fj,s time resolution. 
Effective exposures were 42.5 ks, 10.6 ks, 10.6 ks, and 33.5 
ks with XIS0, XIS1, XIS3, and the HXD, respectively 

Event screening criteria of the HXD are the same as 
Paper I except for two conditions. We allow all geomagnetic 
cutoff rigidity (which was ^6 GV in Paper I), and any data 
transfer criteria, in order to include as many short bursts as 
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possible. These changes do not affect signal to noise ratios of 
short bursts, since they are almost free from the background 
due to their short duration (<1 s). 

We utilized standard filtered files for X1S1 and 
XIS3, while started our XISO analyses from the unfil- 
tered event file. These raw P-sum events were corrected 
and filtered in a way as recommended by the XIS team 
l|Matsuta et al.l I2010T ). including corrections of PI val- 
ues by a ftool software xispt, selections with criteria of 
(GRADE==0GRADE==1GRADE==2) && (STATUS >=0 && STATUS 
and filtering with scripts of xisrepro.xco and xisjrnkf.sel. 
We further eliminated apparent hot pixels from the XISO 
data, and corrected the P-sum mode timing for the known 
delay (31.2 ms in the present case) depending on the source 
location on the CCD chip. 

Figure [2] shows light curves of IE 1547.0-5408 recorded 
with XISO (2-10 keV), XIS1+XIS3 (2-10 keV), HXD-PIN 
(10-70 keV), and HXD-GSO (50-150 keV), which exhibit 
background-inclusive average count rates of 3.89, 5.77, 0.76, 
and 13.38 counts s _1 , respectively. Several short burst events 
clearly appeared in all light curves. 

2.2 Identification of Short Bursts 

The short bursts are characterized as bunched X-ray pho- 
tons. In order to identify such bunched events, we hence 
studied distributions of delta-time (At), which is defined as 
waiting times of individual X-ray photons (including back- 
ground events) from the preceding events. This study uti- 
lized the 10-70 keV HXD-PIN data, since they have a bet- 
ter signal to noise ratio than those of HXD-GSO and XISO, 
and the finer time-resolution (61 /is) than XISO (~7.8 ms). 
Besides, unlike the XIS1 or XIS3 data with burst options, 
HXD-PIN is free from periodic data gaps. The At analy- 
sis provides a standard method to identify bunched events, 
and is actually impleme nted in the HXD digit al electronics 
(HXD-DE; section 4.3 of lTakahashi et alj|2007h to eliminate 
instrumental event bunching on time scales of < a few ms. 

Figure shows the calculated At distribution of the 
10-70 keV HXD-PIN data with a bin size of t bin = 7.8 ms. 
A large portion of these events are non X-ray backgrounds 
(NXB) . The At values of such random events follow an expo- 
nential probability distribution, P(At) oc exp(— XAt), where 
A is the average event rate. The best-fit exponential function, 
shown therein by a dash-dotted line, gives A = 0.778 ±0.001 
count s , which agrees well with the average count rate, 
0.76 count s _1 , derived in t)2.1l On closer inspection of Fig- 
ure |3I a2), the actual At distribution deviates from the best 
fit model below ~50 ms. This indicates that some events are 
bunched, with separations much shorter than A -1 . Thus, we 
choose a threshold interval as Ai t h = 6 tbin = 46.9 ms to 
regard an event as possibly bunched. The chance probability 
of observing At ^ At tn is {1 — exp( — AAt tn )}=3.6%. 

Under the above preparation, we define "an event train 
of length iVbst" as an event bunching where A^st consecutive 
events are detected all with a waiting time of At < Atth- 
Most NXB events sparsely occur with At > At tn , and hence, 
with jVbst=0. Since the chance occurrence of a bunching with 
length JVbst is simply (0.036) JVbBt , we estimate those with 
iVbst =3, 4, and 5 as 4.7 x 10" 5 , 1.7 x 10" 6 , and 6.0 x 10" 8 , 
respectively. Considering the 10-70 keV total PIN events of 
~ 3.4 x 10 4 counts, chance detections of event trains with 
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Figure 4. Distributions of 1-sec count rates for the XIS1+3 
(panel a), HXD-PIN (panel b), and HXD-GSO (panel c) data. 
Best-fit Poisson distributions are also shown. 



Nbet —3, 4, and 5 during the entire observation become ~2, 
~0.06, and ~0.002, respectively. Therefore, event trains with 
Nbst ?S 4 would be securely identified as short bursts. Com- 
pared to a more conventional way of detecting sudden in- 
creases in the light curve, this method is more suited to our 
search for weak bursts on two points; it is unaffected by 
count binning, and it works even when a burt has rather 
gradual rise (e.g., with a precursor). 

Figure |3Jd shows the A^ s t distribution in comparison 
with the predicted chance coincidence probability of ordi- 
nary events, (0.036)^* . We find data points with AT bst > 4 
which are significantly deviated from the fit, and hence re- 
gard them as short bursts. We further merged two separate 
bursts into one if they were detected within the adjacent 1 
sec. Using these criteria, we successfully detected 18 short 
burst events as listed in Table [1] 

To confirm these detections in an independent way, we 
also studied distributions of the count rates in the light 
curves of Figure [5] The calculated histograms are shown in 
Figure [4] In Figure |4jb), the short bursts appear as signifi- 
cant deviations from the best-fit Poisson distribution, which 
is superposed as a solid curve. Even only from Figure [UJb), 
fifteen among the above eighteen bursts were recognized as 
those 1-sec bins where ^8 events were detected. Thus, both 
methods give almost the same detections. Further correc- 
tions for the XIS1 timing mode are described in Appendix IA1 

We also searched the Swift/BAT data for possible detec- 
tions of th e same short bu rsts. Using the same trigger code 
as used in iGrazianil ||2003h on the 64 msec BAT data, only 
2 bursts were recorded by the BAT. This is mainly because 
the BAT, with its wide field of view, has higher background. 
In other words, HXD-PIN onboard Suzaku has a far higher 
sensitivity to weak short bursts, as long as the source is in- 
side its tightly collimated field of view. Therefore, we utilize 
the above Suzaku data set for the following analyses. 
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Figure 5. Continued. 
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3 ANALYSIS 



3.1 Detected Bursts 



Figure [5j shows the XIS, HXD-PIN, and HXD-GSO light 
curves of the bursts described in H2.2I Table Q] also gives 
their basic properties. The bursts are clearly detected by 
HXD-PIN and XISO, and a considerable fraction of them by 
HXD-GSO as well. Due to the burst option, the coverage of 
these bursts with XIS1 and XIS3 (not shown in Figure[5j are 
highly incomplete. This however does not hamper our study, 
because we analyze below the XISO and HXD data of these 
bursts. We exclude Burst-6, presumably the strongest one, 
since a part of the HXD data was lost due to saturation of 
the data transfer. Likewise, Burst-10 is discarded, because it 
fell outside the good time interval (GTI) of the XIS. Thus, 
we hereafter analyze the remaining 16 short bursts. They 
are statistically significant with chance probabilities < 1CP 5 
when combining the 1-sec data of XIS, HXD-PIN, and HXD- 
GSO together (Appendix. [B] 1-sec rates are also shown in 
Table [TJ, and are free from data losses due to dead-time or 
pile-ups (Appendix [Cll . 

The burst peaks in Table[T]are defined as the time of the 
maximum count rate of HXD-PIN, which have higher signal- 
to-noise ratios than the others. As can be seen in Figure 
these peaks can be determined with a typical uncertainty of 
a few bins, or <50 ms. Taking a closer look at individual light 
curves in Figure [SJ some burst photons apparently appear 
outside the duration of iVbst (Table [1} defined in t)2.2l Thus, 
using the XISO (2-10 keV) and HXD-PIN (10-70 keV) data 
co-added together, we manually defined tentative burst du- 
rations to cover whole possible burst X-ray photons. Within 
the above durations, we further calculated the time period, 
T90, starting when 5% of the XISO+PIN photons has been 
detected and ending when 95% has been observed. As shown 
in Figure[6](and also Tablefl]), the T90 values range from ~47 
ms to 1.9 s, with the average of 321 ms. This value is close 
to the mean values derived in recent works of the same ob- 
j ect, T90 =258 ms (|van der Horst et~ai1l2012h and 305 ms 



l|Scholz fc KaspillibT 



where the disagreement may have 
arisen because these works utilized different sets of short 
bursts, with different fluences and different energy ranges. 

Figure [7] represents distributions of burst photon counts 
iV T 90 detected during T 90 . Thus, iV T 90 of the 10-70 keV 
PIN data are distributed from ~5 to ~130. Although At9o 
is sometimes slightly different from iVbst by up to ~50%, the 
difference is negligible in the following analyses. 

For spectral analyses of the individual bursts, the XISO 
and the HXD events accumulated over Tgo were utilized as 
source spectra. The background of XISO was extracted using 
actual events on the CCD chip far away from the source po- 
sition. The HXD-PIN background was not subtracted, since 
the average background rate (~0.8 Hz in 10-70 keV) im- 
plies at most ~1 photon in T90. On the other hand, for the 
HXD-GSO analyses, we produc ed the background sp ectra 
from the simulated GSO events |Fukazawa et alj|2009h dur- 
ing a 20 sec period around the corresponding burst, and sub- 
tracted them after scaling to the duration of T90. The GSO 
background typically contains ~8 counts per bursts, which 
amount to ~30% of the signal counts. Although we did not 
produce background from the actual GSO data taken be- 
fore/after the bursts to avoid possible contamination of the 
persistent or burst emission from the source, this alternative 




Figure 6. Distribution of Tgo for the 16 bursts determined from 
the 10-70 keV HXD-PIN data. 




Log(Counts) 

Figure 7. (a) Distribution of burst counts recorded over T90 
in the 10-70 keV HXD-PIN data, (b) Same as panel (a), but 
for the total counts summed over the 2-10 keV XIS, 10-70 keV 
HXD-PIN, and 50-150 keV HXD-GSO data. 



method gives consistent results. We employed standard re- 
sponse files (epoch 5) for the HXD, while used the response 
(rmf) and auxiliary response (arf) files for XISO produced 
as in 

The average fluxes and fluences of these bursts were es- 
timated using a single PL model, with a column density of 
the photo-absorption fixed at Ah = 3.2 x 10 22 cm~ 2 after 
Paper I. The resultant photon indices and the 10-70 keV 
X-ray fluxes are listed in Table [1] In some bursts, this PL 
model did not give an acceptable fit mainly due to a high 
energy signal deficiency. For such bursts, we alternatively 
utilized a cutoff power-law model (CutPL) with the same 
fixed photo-absorption. The 10-70 keV fluences integrated 
over T90 are also shown in Table [T] As one of the character- 
istics of the present burst sample, Figure [8] shows harness 
ratios of the bursts photons detected by the HXD to those 
by the XIS, as a function of their fluences. The fluences are 
distributed in the range of ~ 10 _9 -10~ 7 erg cm -2 , and their 
hardness ratios ranges from ~0.4 to ~3. Although a hard- 
ening trend toward s weaker bursts was reported from SGRs 
dGogiis et al.ll200ll). and a softening trend was reported from 
AXPs (|Gavriil et all 120041 ), the present sample in Figure [8] 
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Figure 9. X-ray spectra of the three brightest bursts in our sample; Burst-7 (panels al-a5), Burst-14 (panels bl-b5), and Burst-16 
(panels cl-c5). Raw count spectra with XIS0, HXD-PIN, and HXD-GSO are shown in panel (1) with the best-fit CutPL model, and 
corresponding residuals are shown in panel (2). Panels (3) are residuals when the data are fitted by the 2BB model. Panels (4) are 
comparisons among the inferred best-fit CutPL (red) and 2BB (green) models. The spectra are shown after multiplied by individual T90 
values. The XIS0 and HXD-GSO background were subtracted as described in i|3.1l 



10 1 



10 1 

>< 

6 

13 
+ 

Q_ 

o 10° 

■i-l 
fU 
CC 

l/l 
1/1 
CD 
£Z 

~o 

fD 
I 

lO" 1 









H 


12 

H 










16 ^ 

' 9 


5 




i 


ihJj 1 




2| 
17 


p 1 15" 
18 





1CT 9 10 s 10" 7 

Fluence (erg cm" 2 ) in the 10-70 keV 



Figure 8. Hardness ratios of the detected X-ray photons Nxgo 
during T90 calculated as (PIN+GSO)/XIS0, as a function of the 
10-70 keV fluence. Numbers in the figure indicate the correspond- 
ing burst ID. 



does not show any significant corrections between the hard- 
ness and fluence. 

Among the above 16 short bursts, we selected, for our 
detailed spectral analyses, three outstandingly bright ones, 
Burst-7, 14, and 16, which all satisfy iVbst = 60 and show 



fluences above 2 x 10 -8 erg s _1 (Figure |8}. Although these 
are the brightest bursts among the present sample (except 
for Burst-6 which was discarded), they are still weaker com- 
pared to most of the bursts used in previous studies, of which 
the fluences are typically ~10 -8 -10 -4 erg cm -2 . Below we 
apply several spectral modelings to these bursts in the ~1- 
300 keV energy range. For all models, the column density 
for photo-absorption is fixed at Nh = 3.2 x 10 22 cm -2 as 
before, unless otherwise noted. The normalization factor of 
XIS0 to the HXD was fixed at 1.08 based on the correction 
as described in ^A2I although allowing it to float has an 
insignificant effect on the following results. 

The single PL model with the fixed JVh gave acceptable 
fits to none of the three bursts, with reduced xt values of 
4.92, 2.83 and 3.89, for Burst-7, 14, and 16, respectively. 
Even if we make ^Vh free, the fits were still unacceptable 
(xt ~2). Discrepancies of these fits originate from model 
excess in the lowest and highest spectral ends. An optically 
thin bremsstrahlung model is not successful either (xt > 2). 
Then, we fitted these spectral data using the CutPL model, 
and obtained acceptable fits as summarized in Table [2] and 
in Figure [5] The fits yielded a photon index of T ~ 0.1- 
0.8 and a cutoff energy of i7 C ut ~27-66 keV. As alternative 
models, we utilized a two blackbody model (2BB). As shown 
in Table [5J it was also successful on the three bursts, and 
yield the lower and higher temperatures of 35 keV and 1820 
keV, respectively. 
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In Figure [9] (a4), (b4) and (c4), the above two best-fit 
spectral models (the CutPL and 2BB) are compared in vF v 
forms. They give relatively similar spectral shapes, except 
in the higher energy range above ~100 keV where the mod- 
els become unconstrained. In both models, spectral peaks 
appear in the HXD-PIN band (~50-100 keV). Employing 
conventionally the CutPL model, absorption-corrected 0.2- 
300 keV fluences of Burst-7, Burst-I4, and Burst-16, are 
obtained as 1 .1 x 10~ 7 , 9.3 x I0~ 8 , and 8.7 x 10 
respectively. The fluence becomes ~4-8xl0~ 8 
calculated in the 10-70 keV band. 

In Figure [TOk , we compare the spectrum of the brightest 
Burst-7 with the persistent X-ray spectrum recorded during 
the same observation (Paper I). While the persistent spec- 
trum is apparently composed of two components, the burst 
spectrum is more curved, without apparent evidence for such 
two-component nature. 

3.2 Weak short bursts 

As shown in Figure |8j the remaining 13 short bursts have 
considerably lower 10-70 keV fluences than the three stud- 
ied above, distributed below 2 x 10 _s erg cm -2 . They have 
poorer statistics, and also tend to show similar harness ra- 
tios around ~1.0. Therefore, we have stacked their spectra 
together for detailed analysis, with accumulated total expo- 
sure of 3.7 sec. In order to justify the stacking procedure, we 
took spectral ratios of each burst to the stacked one, to find 
that the ratios can be fitted in each case successfully by a 
constant with a reduced chi-square of <1.0. Therefore, the 
13 bursts are concluded to have consistent spectral shapes, 
and hence the stacking procedure can be justified. The de- 
rived 13 constant ratios are distributed from 0.55 to 3.3, with 
the average and standard deviation of 1.56 and 1.05, respec- 
tively. This distribution, ranging by a factor of 6, agrees with 
that of the fiuecne shown in Figure [5] The analysis here uti- 
lized the same responses as those of the brighter three short 
bursts, and in the same way as the previous analyses; the 
background was subtracted from the XIS0 and HXD-GSO 
data. 

Figure [TT] shows a stacked light curve of these weaker 
short bursts accumulated with reference to their peak times. 
Thus, the burst emission is highly significant even in the 
HXD-GSO band. Figure [T2b shows the raw spectrum of this 
cumulative weak-burst data after the background subtrac- 
tion. Its average 10-70 keV flux is by an order of magnitude 
lower than those of the three brightest ones. The HXD-GSO 
background becomes comparable to the signal level around 
~130 keV, and we can claim the HXD-GSO detection at 
least up to 150 keV at 2.8a. 

As summarized in Table a PL model with the fixed 
JVh failed to give an acceptable fit (xt ~ 2.1; Figure [I2p), 
while a PL with free Ah was more successful (xt ~ 1-3), 
yielding T = 1.57 ± 0.04 and A/h = 5.4±°; 8 x 10 22 cm" 2 
(Figure \T2b). In order to further improve the fit especially 
in higher energy range, we again tried the CutPL and 2BB 
fits with the same column density fixed at A^h = 3.2 x 10 22 
cnr 2 . As summarized in Tableland shown in Figure [T2"h , 
the CutPL model gave an acceptable fit (xt ~ 0.8) with 
T = 1.03 ± 0.07 and E cut = 62.9t^;| keV, implying a mild 
spectral curvature. Since V and E cut couple with each other, 
we show in Fig. [13] the fit confidence contours on the F 
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Figure 10. (a) Comparison of X-ray spectra of IE 1547.0—5408 
in vF v forms. Top (black), middle (red), and bottom (blue) spec- 
tra represent Burst-7, the accumulated weak short bursts, and 
the persistent emission from Paper I, respectively. The CutPL 
model was employed to deconvolve the Burst-7 (Table [2]l and 
the accumulated weak- burst spectra (Table |2}, while a PL plus a 
BB for the persistent emission. The column density of the photo- 
absorption was fixed at 3.2 X 10 22 cm -2 in all cases, (b) The ratio 
between the red and blue spectra in panel (a), (c) Same as panel 
(b), but after eliminating the BB component from the persistent 
emission. 



vs. -Ecut plane. In contrast to the successful CutPL model, 
the more convex 2BB model, which was successful on the 
brightest three bursts ( H3.ip . became much less successful 
(xt ~ 1-6; Figure [12b). Thus, the weaker bursts are consid- 
ered to have a flatter 10-70 keV HXD-PIN spectrum than 
the brightest bursts, particularly Burst 7 and 16. 

To make the above spectral difference clearer, we added 
this stacked weak bursts to Figure [lOk in vF v form, where 
we employed the CutPL model for deconvolution in the same 
way as Burst-7. The weaker bursts are by ~2 orders of mag- 
nitude brighter than the persistent emission, and by ~1 
order of magnitude fainter than Burst-7. As visualized by 
this plot, the cumulative burst shows a hard X-ray spec- 
trum which is less curved than that of Burst 7 and is sim- 
ilar to that of the persistent X-rays. In fact, the value of 
Tbst = 1-57 ± 0.04, obtained above by the PL fit with free 



Ah, is consistent with r per = 1.54 



+0 03 of the persistent 



hard component (Table [3]). Although the CutPL model gave 
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Figure 11. Light curves obtained by stacking the 13 weaker 
short bursts, with XIS0 (panel a), HXD-PIN (panel b), and HXD- 
GSO (panel c), in the 2-10, 10-70, and 50-150 keV band, respec- 
tively. 



a much harder photon index, Tbst = 1-03 ± 0.07, the slope 
is still similar, if considering the effect of cutoff around ~63 
keV. 

In order to more directly compare these vF v spectra, 
we divided the stacked weak-burst spectrum by that of the 
persistent emission. As shown in Figure 110b . the resultant 
ratios stay constant at ~170 over the ~8-200 keV energy 
range. The remaining difference in the ratio below ~8 keV in 
Figure [TOb is thought to originate mainly from the presence 
of the soft BB component. We hence re-calculated the same 
ratio in Figure UOfc after eliminating the soft BB component 
from the persistent spectrum. Thus, the ratio became much 
flatter even in energies below ~10 keV, suggesting that the 
weak-burst spectrum is similar in shape to the hard X-ray 
component in the persistent emission. 

Finally, in order to examine our weak-burst sample for 
its homogeneity, we subdivided the 13 bursts rather ran- 
domly into two subsets, under a constraint that they should 
be comparable, to within 10% in the summed signal photon 
counts (XIS0, PIN, GSO summed). From the two subsets, 
we then produced two stacked spectra instead of one, and 
analyzed them in the same manner as before. The pair of 
CutPL parameters obtained in this way are indicated in Fig- 
ure [13] by a pair of data ponits with the same symbol. By 
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Figure 12. (a) Background-subtracted and response-inclusive 
spectra of the stacked weak bursts from IE 1547.0—5408, com- 
pared with the best fit CutPL model (Table [3J. Panels (b)-(f) 
show residuals when using a PL model with the fixed column 
density, a PL with the column density left free, a CutPL model, 
a 2BB model, and a 2BB+PL model, respectively. 



trying six different partitions, we confirmed that the two 
subsets always give, within errors, consistent model param- 
eters, which in turn are also consistent with those derived 
from the entire sample (indicated by contours). Therefore, 
the present sample can be regarded, within errors, as homo- 
geneous. 



3.3 Fitting by thermal/non-thermal components 

Let us further examine the suggestion of Figure [10] that the 
burst emission has the same spectral components as the per- 
sistent X-ray emission. Especially the good similarity above 
~8 keV points to possible presence of the hard PL compo- 
nent in the stacked weak-burst spectrum, like in the persis- 
tent emission spectra. The short bursts from the transient 
magnetar SGR 0501+4516 were already e xamined for such 
a po ssibility with an affirmative answer (|Nakagawa et al] 

Eon]). 

In order to assess the above consideration, we re- 
analyzed the stacked weak-burst spectrum, by adding, to the 
2BB model, a hard PL with its slope fixed at r pcr = 1.54 
as specified by the persistent emission (Paper I). The inclu- 
sion of the hard PL has made the 2BB fit acceptable, as 
summarized in Figure I12F and Table [3] The additional PL 
contributes to the total 1-300 keV flux by 52%, and is sta- 
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Figure 13. Confidence contours at the 68% (light blue), 90% 
(green), and 99% (red) confidence levels of the CutPL (panel a) 
fit to the stacked weak burst spectrum. Individual data points 
show the parameters, with ±lcr errors, when the 13 short bursts 
are divided into two groups in six different ways. 



tistically significant, because the achieved fit improvement, 
A\ 2 = —27.1, implies an F-test chance probability of of 
2.4 x 10 -7 . Even though an acceptable fit with a compa- 
rable goodness had already been obtained with the CutPL 
model, we consider the 2BB+PL modeling more appropri- 
ate, for the following two reasons. One is that this modeling 
has successfully explai ned the spectra of wea k short bursts 
from SGR 0501+4516 (|Nakagawa et al.ll201ll ). The other is 
that the 2BB model has been regarded as a standard repre- 
sent ation of short-burst spec tra of a fair number of magne- 
tars ([Nakagawa ct al. 20 07bl and references therein), so that 
our final model should also be based on it, rather than the 
CutPL model which in the present case is quite conventional. 

Now that evidence for the PL component has been ob- 
tained in the stacked weak-burst spectrum, it would be a 
natural step to examine the brighter ones, Burst-7, 14, and 
16, for the same possibility, even though their spectra were 
already represented successfully by the 2BB model. This 
is because the distinction between the brighter and fainter 
bursts in the present observation is purely technical, with- 
out particular reason to consider that the two groups have 
distinct nature. Then, to the 2BB model describing the 3 
bright bursts, we added the hard PL with Pbat fixed at the 
same value. This improved the fit to Burst-7, 14, and 16, 
by Ax = —4.3, —2.0, and —2.8, giving F-test probabilities 
of 0.60%, 7.4%, and 2.4%, respectively. Thus, the hard PL 
component is likely to contribute at least to Burst-7. The 
best fit spectra are shown in Figure [T3] in vF v forms, with 
the derived parameters in Tabel [2] The BB temperatures 
did not change within errors from those obtained in the pre- 
vious pure 2BB modeling ( H3.1[l . while the hard PL compo- 
nent was found to contribute to the total 1-300 keV fluxes 
of Burst-7, 14, and 16, by 10%, 14%, and 19% respectively. 



(a) Burst-7 



(b) Burst- 14 





(d) Weak Bursts 



Figure 14. An vF u presentation of the 2BB plus PL (with 
r = 1.54 fixed) fits to the spectra of Burst-7 (panel a), Burst-14 
(panel b), Burst-16 (panel c), and the stacked weak short bursts 
(panel d). The two BB components are represented in blue and 
green, while the hard PL in red. 



4 DISCUSSION 

4.1 Spectral comparison among different 
emissions 

During the 33.5 ks of effective exposure with the HXD onto 
IE 1547.0-5408 in its 2009 January activity, we detected 18 
short bursts with their fluence S above S m i n ~ 2 x 10 -9 erg 
cm -2 , and analyzed 16 of them. Even the strongest three 



among the 16 bursts, with the 10-70 keV fluence of S = 4- 
8 x 10~ 8 erg cm~ 2 (or ~ 1 x 10~ 7 erg cm -2 in 1-300 keV), 
are still weaker than most of previously reported bursts from 
magnetars in general which typically have S > 10~ 7 erg 
cm -2 . For example, the burst forest of IE 1547.0-5408 ob- 
served on 2009 January 22 with Swift, INTEGRAL, Suzaku, 
and Fermi was composed of event s with S > 2 x 10~ 6 er g 
s _1 up to 2.5 x 10~ 4 erg cm~ 2 (|Mereghetti et all I2009T ). 
Therefore, the present Suzaku results, together with those by 
INakagawa et all (|201ll ) on SGR 0501+4516, provide valu- 
able information on the wide-band spectra of weak short 
bursts. For reference, the present sample defines a burst fre- 
quency of 5.4 x 10 -4 burst s _1 in the range of S > S'min- 

So far, wide-band spectra of many energetic (S > 
10~ 7 erg s _1 ) short bursts from magnetars have been 
explained successfully by the 2BB model: these include 
numerous bursts from SGR 1806-20 and SGR 1900+14 
(|Nakagawa et alj|2007bl ; llsrael et~ai1 |200S1) . as well as the 
strongest Suzaku burst from SGR 0501+4516 (|Enoto et al.l 
120091 ; INakagawa et al.l l201ll ). At the same time, the 2BB 
model can generally reproduce the soft X-ray component, 
which is ubiquitously seen at energies of < 10 keV of 
wide- band persistent spectra from magnetars (|Mereghettil 
|200S|) including IE 1547.0-5408 itself (Paper I). These 
2BB fits to the two different forms of magnetar emis- 
sion reveal an interesting co mmon scaling as TLow/Imgh ~ 
0.4 (jNakagawa et all |2009| ). where T Lo w and T Hig h de- 
note the lower and higher 2BB temperatures, respectively. 
These results reinforce the possibility that the persis- 
tent emission is composed of numerous micro short bursts 
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(CTlThompson fc Duncanll 19961 ; iLvutikovllioO^ ; iNakagawa I 
l2007al ). Since the 2BB modeling is thus considered to have 
some physical meanings beyond a mere convention, below we 
construct our discussion based on the 2BB modeling, rather 
than those with the CutPL model which is rather empirical. 

The three brightest bursts in our present sample, par- 
ticularly the brighter two (Burst 7 and 16), among our sam- 
ple exhibit upward convex spectral shapes, which have been 
successfully described by the 2BB model over an extremely 
wide energy band of ~ 1 to ~ 300 keV ( 93.1J I. Importantly, 
the values of TLow and Tmgh derived from these bursts (pos- 
sibly except for Burst- 14) obey the above scaling relation. In 
contrast, the stacked weaker bursts have a less convex shape 
(Figure [TO)) , and did not accept the simple 2BB modeling. 
This is understandable if we presume that the spectra of 
weaker bursts are contributed by the hard component with 
a power-law photon index F ~ 1, which i s generally observed 
in the persistent em ission of magnetars (|Kuiper et al.|[200rj ; 
lEnoto et al1l2010d '). The weak bursts from SGR 0501+4516, 
when stacked together, actually exhibited a hard-tail feature 
l|Nakagawa et all 1201 ll ). In fact, Fi gure [10] b and c provide 
supporting evidence for this possibility, where the ratios be- 
come flatter above ~8 keV. We have accordingly added a 
hard PL component to the 2BB model in A3. 31 fixing V at 
the same value as obtained from the persistent emission, 
and obtained a fully acceptable spectral fit. Thus, the hard 
component is significantly present in the stacked weak-burst 
spectrum, as long as the 2BB model is chosen as the start 
point of analysis. As a further justification of this 2BB+PL 
modeling, the stacked spectrum yield Tmgh/TLow = 4.9±1.1 
(from Table[3| in agreement with the "canonical" ratio, with 
it was 7.3 + 0.5 before properly i nclud ing the PL component. 
As argued bv INakagawa et all ()201ll ). these strengthen the 
similarity between the burst and persistent emissions. The 
report by llsrael et al.l |2008l ), that weaker bursts tend to 
show harder spectra in terms of the 2BB analysis, may be 
understood if considering possible contributions of the same 
PL component to their weaker bursts. 

The above results have lead us to a natural question 
about brighter short bursts, of which the spectra generally 
lack the PL component: is this an intrinsic effect due, e.g., 
to some changes in the emission mechanism, or an artifact 
caused by their higher 2BB temperatures which mask any 
PL emission? With this in mind, we examined ( t]3.3[) the 
spectra of the three brighter Suzaku bursts, for possible con- 
tribution from such a hard component. Then, the presence 
of the PL component (on top of the 2BB model) was sug- 
gested in Burst-7 at a relatively high confidence (>97%) 
level. The results are presented in Figure [15] on the ther- 
mal (2BB) vs non-thermal (PL) luminosity plane, in com- 
parison with the weak-burst and persi stent spectra. There , 
assuming a source distance of 4 kpc ()Tiengo et al.ll2010h . 
the persistent X-ray emission (Paper I) is represented by a 
bolometric blackbody luminosity of Lbb = 1-2 x 10 35 erg 
s _1 and absorption-corrected 1-300 keV PL luminosity of 
Lpl = 6.9 x 10 35 erg s , which sum up to give the to- 
tal luminosity of L peT = 8.1 x 10 35 erg s - 1 . We also quote 
the re sults on SGR 0501+4516, taken from INakagawa et al.l 
l|201ll) . Although it is still difficult at present to distinguish 
the two alternatives the PL luminosity, if it is universally 
present in short bursts, appears to saturate at ~ 10 38 erg 
s _1 , preferring the former possibility. 
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Figure 15. Comparison of the 1—300 keV PL luminosity to the 
bolometric BB lumin osity. The results on SG R 0501+4516 were 
taken from Figure 8 in Nakagawa et al. (201l|). In this plot, errors 
are shown at the 90% C.L. level. 



On large scales, Figure [TS] reveals a clear positive corre- 
lation between luminosities of the two emission components. 
Therefore, the soft (thermal) and hard (PL) components are 
considered to approximately keep their luminos ity ratio over 
a very broad scale in flux or fluence. Similarly, lEnoto et al.l 
(|2010ch showed that the luminosity ratio between these two 
emission components, comprising the persistent emission, is 
not much different between magnetars in activity and those 
in quiescence. These facts give another support to the mi- 
croburst conjecture. 



4.2 LogN-logS relation of short burst 

In order to examine the possibility that the persistent X-ray 
emission is composed of micro burst events, it is inevitable 
to evaluate X-ray fluxes accumulated over resolved and un- 
resolved bursts events. Since the present short bursts, with 
S > Smin ~ 2 x 10 -9 erg cm~ 2 , were detected at a low 



at a level of ~ 10~ erg s~ cm~ . This is by ~2 orders 
of magnitude lower than the observed persistent X-ray flux. 
Therefore, much larger contributions of smaller unresolved 
short bursts are required to explain the persistent X-ray flux. 

For a quantitative estimate, let N(> S) denote the oc- 
currence frequency of those bursts of which the fluence is 
> S. Observationally, N is thought to be described by a 
single power-law as N(>S) oc S~ a , like in solar flares and 
earthquakes. Sometimes known as the Gutenberg-Richter 
law, this relation is considered to represent self-organized 
criticality. Previous observations of SGR 1806—20 gave a 
range of a = 0.7-1. 1; e.g., a ~ 0.9 by the Konus- Wind 
|Aptekar et al.l 1200 ll ). a = 0.76 ± 0-17 by the BATSE, 
a = 0.67 + 0.15 b y ICE dG5gus et al.ll2000h , a = 0.91 + 0.09 
by INTE GRAL ilGotz et 8.1.1120061 ), and a = 1.1 ± 0.6 by 
HETE-2 (|Nakagawa et al.l l2007bf ). From the present tar- 
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Figure 16. Expected luminosity of integrated unresolved weaker 
bursts, as a function of the size-distribution slope a. Blue, green, 
and red lines are the calculations with the lower bounds of the 
integration at Si/S m i n = 10~ 2 , 10 -3 , and 10~ 5 , respectively. 
Three horizontal lines represents the persistent X-ray luminos- 
ity of IE 1547.0-5408 for the total (top), PL (middle), and BB 
(bottom) components. 



get, IE 1547.0—5408, the slope is reported to be a = 
0.75 ± 0.06 and a = 6 ±0.1 through I NTEGRAL/ACS 
jMereehetti et all 120091 ) and Swift/XRT (jScholz fc Kaspil 
l201lj ) observations, respectively. 

The burst number density at a particular fluence S is 
written as f hst (S) oc -dN(S)/dS = (f /S a )(S/S )- a - 1 , 
where Sq is a typical fluence and /n (burst s -1 ) is a nor- 
malization constant. Thus, the luminosity integrated over 
bursts between Si and 52, at a distance d (kpc), is given as 



S-2 



4nd 2 Sf bat (S)d(S) 



4nd J Sofo 



4-Kd 2 S f (si, 



a ds 



s{- a )/(l-a), 



Si/So, 



and S2 

9 



S2/S0. Here, let 
2 (the present 



with s = 5/5*0, Si 

us assume So = S m i n ~ 2 x 10" " erg cm 
detection limit), and S2 = S2/S m i n = 1, at d =4 kpc. The 
observed frequency N(Smin) ~ 5.4 x 10 -4 burst s _1 gives 
/o ~ 7 x 10~ 4 burst s _1 . In order to explain the persistent 
luminosity in this way, we clearly need a > 1.0, because 
S x N oc S 1 '" should increase towards smaller values of S. 
More quantitatively, Figure \W\ shows the calculated value of 
Lx as a function of a for these different values of si. 

Comparing the three brightest bursts and the stacked 
weaker ones in Table [2] and Table [3] the blackbody radii 
of neither the lower- or higher-temperature components are 
strongly dependent on their fluences S. Thus, as a first-order 
approximation, let us assume S oc kT 4 after the Stefan- 



Boltzmann law. Since the weak bursts with S m i n and persis- 
tent emission give fcTkigh ~ 12 keV and kT ~ 0.65 keV, re- 
spectively, a typical unresolved weakest burst is considered 
to have a fluence si = Si/S min ~ (0.65 keV/12 keV) 4 ~ 
10~ 5 . Combining with Figure [TBI this implies that the slope 
a should be ~1.2 and ~1.4 to supply the luminosity of the 
persistent BB and PL luminosities, respectively. 

The required value of a = 1.2-1.4 is larger than has 
been measured. Since the observations of bursts were per- 
formed mainly during more burst-active states than in the 
Suzaku observation, a possible interpretation is variations 
of a depending on the magnetar activity. For example, the 
solar flare is known to change its size-distribution slope a, 
from ~0.8 at act i vity p eaks to ~1.2 in more quiescent states 
l|Bromund et al.l Il995l ). This is also pointed out in some 
magnetars (see also S CR 1806-20 and SCR 19004-14 in 
iNakagawa et alj2007bl, r elated with SCR 0501+4516 behav- 
ior in Gogtis et a.1.1 l2O10h ■ If a reaches 1.2-1.4 in less burst- 
active states, the persistent emission can be explained by 
the superposition of micro bursts. Otherwise, an alternative 
dissipation process that is independent of the short-burst 
production is inferred to be also contributing to the persis- 
tent emission. Such p rocesses include, e.g. , internal crustal 
heating in magnetars l|Aguilera et aLlfeOQct ). 



4.3 Energy source for the emissions 

Finally, let us compare energetics of the bursts and per- 
sistent emission with the magnetic energies of magnetars. 
Assuming 4 kpc as before, the 2BB model and the Stefan- 
Boltzmann law measure the total radiated energies of 2.1 x 
10 38 erg, 1.6 x 10 38 erg, and 1.6 x 10 38 erg, from Burst-7, 
Burst-14, and Burst-16, respectively. Further considering the 
Tgo values, these yield bolometric luminosities of 7.8 x 10 38 
erg s~ 4 , 1.9 x 10 38 erg s -1 , and 5.9 x 10 38 erg s _1 , respec- 
tively. In the weaker bursts, an average energy emitted per 
bust becomes 1.2 x 10 37 erg, and the bolometric luminosity 
is 4.2 x 10 37 erg s~\ 

In comparison with the above estimates, the magnetic 
energy stored within a size of -Rmag becomes 



-Ebst 



^mag *^g 

2 x 10 



B 



B s 



-Rns 



erg, 



when norma lizing to the dipo le magnetic field, -B sur f = 
2.2 x 10 14 G (|Camilo et al.ll2007t ). and employing the canon- 
ical neutron start radius _Rns = 10 km. If the individual 
burst emission is powered by dissipation of the magnetic en - 
ergy in this region (e.g., via reconnection; iLvutikovl l2006t ) ■ 
a typical size of i? mag > 10 _2 ' 3 7?ns is needed. Blackbody 
radii observed in our sample (the 2BB model), i?Low~l km 
= 0.1-Rns and i?m g h~0.1 km = 10~ -Rns, are larger than 
the required i? mag . This is considered reasonable, because 
the released magnetic energy would in any case diffuse out, 
and would make Rbb > flmag- 

If we assume si ~ 10~ 5 and a ~ 1.4 like in the previ- 
ous subsection, the burst frequency at Si is estimated to be 
/ ~ 10 3 Hz. Therefore, a minimum size of -R ma g /Rns ^ 10 
is enough to supply the observed luminosity L x ~ 6 x 10 35 
erg s _1 , if the released magnetic energy is evaluated as 
(B/B sult ) 2 (R mag /R NS ) 3 f erg s" 1 . Since the 
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persistent emission exhibits -Rbb ~ 5 km (Paper I) , the con- 
dition of Rbb > -Rmag is retained. 



5 CONCLUSION 

We studied short bursts from IE 1547.0—5408 (also known 
as SGR J1550-5418, PSR J1550-5418) detected during the 
Suzaku pointing observation on 2009 January 28-29 (UT) in 
its burst-active state. This period was 7 days after the burst 
forests on January 22. Combined with the previous study 
of its persistent X-ray emission (Paper I) , we have obtained 
the following results. 

(i) Using a At distribution of the HXD-PIN data, we 
identified 18 short burst events with its risk probability be- 
low ~ 1.7 x 10 -6 . Sixteen of them, free from the data loss, 
were used for the spectral analyses. 

(ii) Three brightest bursts have their fluences of ~ 4- 
8 x 10~ 8 erg cm~ 2 in the 10-70 keV energy range. Their 
individual spectra were fitted successfully by the 2BB model 
(and some other models) over the 0.5-^400 keV band. 

(iii) Remaining 13 bursts define one of the weakest sam- 
ples ever measured, with their fluences covering the range of 
2 x 10 _9 -2 x 10 -8 erg cm -2 . Their stacked spectrum shows 
similarity to the persistent emission in its slope above ~8 
keV. This spectral similarity extends to even towards lower 
energies after eliminating the BB component from the per- 
sistent spectrum. 

(iv) The stacked spectrum cannot be represented by a 
2BB model, while an additional hard PL made the fit ac- 
ceptable, even fixing the slope at the value of the persistent 
one. The luminosity of the hard PL shows a sign to saturate 
at ~10 38 erg s _1 , when compared to the thermal luminosity. 

(v) We evaluated the available energy supplied from 
the unresolved short bursts. The persistent emission from 
IE 1547.0—5408 can be explained if the slope a of the cu- 
mulative distribution of the short bursts is rather steep as 
a = 1.2-1.4. 

We thank members of the Suzaku magnetar Key Project 
and the Suzaku operation teams for the successful ToO. TE 
is supported by the JSPS Postdoctoral Fellowships for Re- 
search Abroad. 
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Table 1. List of detected short burst events. 



ID 


Time at Burst Peak 


T90 


N L 


N 2 


1-sec Rate 2 


Model 


Photon 


Flux 3 


Flucncc 3 






(UTC) 


(ms) 






Index 






(p- value) 


1 


2009-01-28T21:48:24.8 


1851.4 


9 


(72, 26, 29) 


(37, 15, 25) 


PL 


1.65+S"" 

— U. lo 
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9.0+1-1 


1.09(0.36) 
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2009-01-29T02:36:02.2 


46.9 
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(16, 9, 4) 


(26, 10, 25) 


PL 


1.34+°'" 


6.0+\ 6 


2.8 +0 5 T 

— U. ( 


0.81(0.60) 


3 


2009-01-29T03:06:21.3 


148.4 


5 


(9, 4, 7) 


(20 , 6, 14) 


PL 


l-°7 + : 25 
— u.zo 


1.3+° : ? 


i.9+;-; 


0.75(0.69) 


i 


2009-01-29T08:09:06.6 


226.5 


4 


(17, 8, 10) 


(21, 10, 21) 


PL 


— U. 12 


1.5+° : 1 

— U.4 


3-4+o : o 
— u.y 


1.19(0.31) 


5 


2009-01-29T08:53:23.4 


78.1 


5 


(8, 5, 4) 


(12, 9, 14) 


PL 


— U.Li 


24 +o.5 


1.9+ - 4 

— 0.8 


0.44(0.82) 


6 


2009-01-29T09:17:10.4 


656.2 


256 


(2029, 219, 59) 


(2107, 282, 93) 


PL 








- 


7 


2009-01-29Tll:19:27.4 


265.6 


139 


(127, 131, 84) 


(155, 143, 101) 


CutPL 


u - lz -0.14 


4-0 7 

28-iiS:I 




0.82(0.74) 


8 


2009-01-29T12:16:56.4 


390.6 


19 


(44, 23, 15) 


(56, 28, 26) 


PL 


1 49+0-08 






1.46(0.15) 


9 


2009-01-29T12:17:36.5 


54.7 


11 


(19, 11, 10) 


(34, 16, 32) 


PL 


QQ1+0 .13 

u - 9± -0.13 




8.2+J;° 


0.28(0.95) 


10 


2009-01-29T12:58:41.1 


78.1 


11 


(-, 38, 29) 


(-, 47, 41) 


PL 










11 


2009-01-29T13:41:57.1 


164.1 


8 


(21, 8, 11) 


(37, 10, 20) 


PL 


1 45+0.14 


n c+0.5 
z -°-0.6 


4 1+0.8 
4ii -0.9 


0.85(0.55) 


12 


2009-01-29T13:43:27.9 


210.9 


7 


(7, 9, 13) 


(13, 10, 26) 


CutPL 


n no+0.39 


4 o+°- 8 
^■"-0.7 


8-411:1 


0.53(0.78) 


13 


2009-01-29T13:46:37.4 


328.1 


13 


(97, 20, 19) 


(113, 27, 26) 


CutPL 


1 li+o- 17 

J - J " L -0.18 


4 3+0-2 
^• J -1.3 


u.o±H 


0.87(0.56) 


11 


2009-01-29T13:50:38.2 


835.9 


60 


(156, 62, 82) 


(187, 67, 89) 


CutPL 


n OO+0.09 


5 rr+0.3 
°-°-0.8 


45.9t 2 -j 


1.10(0.34) 


15 


2009-01-29T15:40:45.1 


109.4 


15 


(35, 14, 14) 


(44, 16, 24) 


PL 




4 8 +0 - 7 


5 O+0.7 


1.23(0.26) 


16 


2009-01-29T16:51:40.5 


265.6 


90 


(97, 85, 75) 


(148, 96, 98) 


CutPL 




21 1+ - 9 


56.ll 2 7 j 


0.54(0.97) 


17 


2009-01-29T16:59:58.8 


109.4 


6 


(19, 6, 5) 


(24, 8, 19) 


PL 


I QO+0.17 
1 - oy -0.15 


2.6+S-J 


2 9+ - 8 
z.y_ 1 2 


0.78(0.60) 


18 


2009-01-29T18:42:12.6 


54.7 


5 


(17, 6, 7) 


(25 , 8, 23) 


PL 


i OO+0.19 




3 6+ - 9 


0.27(0.97) 



1 : An event train of length JVbst is defined in i]2.2l 

2: iVx90 arl d 1-sec rate around the burst are shown in a form of (XIS0, PIN, GSO). The former is a burst photon count during Tgo 
defined in i|3.1l The latter is extracted from 1-sec binned light curves in Figure [2] 

3: Flux is shown in a unit of 10 — 8 erg s _1 cm -2 , and flucncc is in 10 — 9 erg cm -2 both in the 10—70 keV band. 
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Table 2. Spectral parameters of the three brightest short bursts. 



Spectral Model 


Parameter 


BST-7 


BST-14 


BST-16 




Time (2011-01-29) 
T90 (sec) 

Fluence 1 (10 — 8 erg cm -2 ) 


11:19:27.4 
0.2656 
7 46+ ' 48 


13:50:38.2 
0.8359 
4 60+ - 25 


16:51:40.5 

0.2656 

5 60+ ' 24 
O.DU_ 72 


wabs 2 


Nu (10 22 cm^ 2 ) 


(3.2 fix) 


(3.2 fix) 


(3.2 fix) 


CutPL 


Photon Index T 

Cutoff Energy E cut (keV) 

PIN Flux^(10" 7 erg s" 1 cm" 2 ) 

Total Flux^lO -7 erg s" 1 cm -2 ) 

fit goodness x 2 {y) 

figure 


12+ ' 13 
u - lz -0.14 

26.9+*-| 

2 -8 1+ 03i 

4 01+ 014 
i.ui_ 51 

0.82 (28) 

figl2tal,2,5) 


0.83 ± 0.09 

0*^^ — 08 

1.10 (20) 
fig|§tbl,2,5) 


0.28 ±0.12 

3L9l 4 ' 8 

n 11+0.09 
1 - ± q 27 

o 91+O.II 

°- zl -0.56 

0.54 (25) 
figEIcl,2,5) 


2BB 


feT Low (keV) 
fiLow (km) 
fcT High (kcV) 
^High (km) 

T3T1\T T?li.-..5/in — 7 nvrr n — 1 ,. .... — 2 "\ 

r^llN r lux ^1U erg s cm ) 
Total Flux^(10~ 7 erg s" 1 era" 2 ) 
fit goodness %v { v ) 
figure 


r- 9 n+0.43 
OiZ °-0.42 
i 7s +0.20 

20.5 ±2.1 
142 +0.038 
-0.028 
9 r-,,+0.12 
2 ' 54 -0.35 

405 +o.i3 

0.62 (27) 
ng|Ha3,5) 


„ to+0.47 
z - oo -0.34 

20.9+J; 7 
0.0766j:°;°?J« 
0.45_ 12 
0.98 ± 0.02 
0.79 (18) 
fig|9lb3,5) 


3.97 ±0.41 

9 1O + 0.31 

2.18_ 27 
17.5+i' 8 

191+0 047 
u - 13 -0.0381 

1 qi+0.12 
^ -0.29 
05+0- 12 
d - u0 -0.45 

0.58 (24) 
figHc3,5) 


2BB + PL 


fcT Low (keV) 
4ow (km) 
kT msh (kcV) 
Rm s h (km) 

Tbst 

PL Flux^lO" 7 erg s" 1 cm" 2 ) 

fit goodness xi (") 

figure 


c or+0.67 
°- oo -0.57 

1 45+ ' 24 

20.2 +2 j 

1QO+0 042 
u. lo»_ 033 

1.54 (fix) 

410+° 190 
u.iiu_ lg5 

0.48 (26) 

figda) 


2 90+ - 30 

z - au -0 .28 
1 QQ+0-37 
±,a -0.33 

7Q.9±\% 

0.0730+ 9 ; , 2 " 

1.54 (fix) 

140+ ' 092 
U.±1U_ 093 

0.68 (18) 

figlHb) 


4 0,5+0.63 

r^+0.050 
U.ODO_ 04g 

17 1+ 1 - 9 
191+0.052 
u.iyi_ 04i 

1.54 (fix) 

ooo+0. 196 
u. oo»_ 200 

0.49(23) 

figLHc) 



1: Fluences are estimated using the CutPL model in the 10-70 keV range with T90. 
2: The interstellar absorption (wabs) is multiplied to each model. 

3, 5, 7: Absorbed c X-ray fluxes in the 10-70 keV energy range at its average value. 

4, 6, 8, 9: Absorbed X-ray fluxes in the 1—300 keV energy range at its average value. 
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Table 3. Spectral parameters of the accumulated weaker short burst. All the quoted errors are at the la level. 



Burst Emission 




Summed exposure (sec) 
Ave. fluence 1 (erg cm - 2 ) 


3.7 

6.61+°;°* x 10 _ 9 


Spectral Model 


Parameter 


Value 


PL (N H fixed) 


Nh (10 22 cm" 2 ) 
Photon Index r bat 
fit goodness xi ( v ) 


3.2 (fix) 
1.5 

2.05 (34) 


PL 


iV H (10 22 cm" 2 ) 
Photon Index T^t 
Flux 1 (10 -8 erg s~ J cm" 2 ) 
fit goodness x 2 {y) 


5 4+0.8 

1.57 ±0.04 
1.7 ±0.1 
1.29 (33) 


CutoffPL 


N H (10 22 cm" 2 ) 
Photon Index r^at 
Cutoff Energy E cut (keV) 
Flux 1 (10~ 8 erg s _1 cm -2 ) 
fit goodness xi iy) 


3.2 (fix) 

1.03 ±0.07 
62 Q+ 14 - 5 

O n + 0.1 

z -°-0.2 

0.81 (33) 


2BB 


iV H (10 22 cm" 2 ) 
kT Lom (keV) 
JtLow (km J 
feT H i gh (keV) 
Rm g h (km) 

Flux 1 (10~ 8 erg s _1 cm -2 ) 
fit goodness x 2 (f ) 


3.2 (fix) 
1 ' b8 -0.1S 

Q QQ"| - 0.52 

O.OO_ 4 g 
19 q+0.9 

140+ ' 021 
u.rtu_ 01g 

2.5 ± 0.2 

1.63 (32) 


2BB + PL 


N H (10 22 cm" 2 ) 
fcT Low (keV) 
.Rlow (km) 
feTnigh (keV) 
Rm g h (km) 

r h 

PL Flux 1 (10" 8 erg s" 1 cm" 2 ) 
fit goodness xi ( y ) 


3.2 (fix) 
9 71 +0.50 
z - ' -0.48 
1 16+ 041 

13.3± 2 ;° 
0.0905+°S 
1.53 (fix) 
0.78 ±0.13 
0.82 (31) 


Persistent Emission 2 


BB + PL 


Nh (10 22 cm" 2 ) 
kT (keV) 

Tpcr 

Flux 1 (10~ 10 erg s" 1 cm" 2 ) 
fit goodness x 2 (f) 


3.2 ± 0.1 

U - D0 _0.01 
I 54+O.O3 

1.08 (278) 



1: Average fluence is evaluated in the 10—70 keV band assuming the 2BB+PL model. Flux are estimated in the 10—70 keV. 
2: Values from Paper I. The quoted errors are converted to the lcr level. 
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APPENDIX A: CORRECTION FOR THE 
P-SUM MODE XISO 

Here we describe further calibrations of the time assignments 
and responses of the P-sum (timing) mode XISO, dedicated 
towards the present analyses. 



R xy (i) = 1 ^— J2 sO'Mi+i) (* = 0,1,»-,JV-1).(A1) 

Figure IA2f a) and (b) show cross-correlation coefficients of 
the 10-14 keV and 40-70 keV data, respectively. As clearly 
shown in these figures, no time lag was confirmed between 
the HXD 40-70 keV data (Figure lA^b). while the cross- 
correlation peak of the 10-14 keV range appears at 31.25 
ms in Figure lA^ a) . Thus, we assigned the time correction 
of 31.25 ms to the XISO data set (TIME-0.0325 s). This cor- 
rection is consistent with the previous analyses of the per- 
sistent X-ray emission at §3.1 in Paper I. Individual light 
curves of XISO in Figure [5] were already corrected on this 
time lag of the 31.25 ms. 



Al Timing correction 

It is known that there is still a slight time lag (~30 ms) as- 
signed to the P-sum mode of the XIS, even after including 
corrections related with read-out time delays , depending on 
the s ource position on the XIS CCD chip fsee lMatsuta et all 
120091 ). To study the burst light curves in detail, we further 
assigned an additional timing correction to the P-sum mode, 
comparing the HXD-PIN data, w hich was already acc urately 
calibrated using the Crab pulsar l|Terada et alj|2008h . In or- 
der to evaluate the residual time lag, we compared the burst 
light curves of the XISO to that of HXD-PIN both in the 10- 
14 keV. We accumulated 5 short bursts in which the 10-14 
keV photons were clearly detected in the HXD-PIN data. 
Panels in Figure IA1I (left) are the 10-14 keV cumulative 
light curves of XISO and HXD-PIN, respectively, while pan- 
els in Figure \M] (right) are the 40-70 keV ones of HXD-PIN 
and HXD-GSO as a reference. As compared here, the 40-70 
keV peaks correspond with each other, while a peak of the 
10-14 keV XISO data comes slightly later than that of the 
HXD-PIN. 

In order to quantitatively evaluate this XISO time 
lag, we calculated cross-correlations of them. The cross- 
correlation coefficience R xy (i) for the i-th bin is calculated 
by 



A2 Response correction 

For the XISO spectral analyses, we produced rmf and ar f 
files based on the procedure stated in lMatsuta et al] (|201Cf ). 
Since the P-sum spectral data is not widely utilized so far, it 
requires a careful calibration, we further included correction 
factors to rmf and arf fiels, as follows. The persistent X-ray 
emission from IE 1547.0—5408 was already measured using 
the XIS1, XIS3, HXD-PIN, and HXD-GSO data in Paper I. 
We also produced the same persistent X-ray spectra from the 
XISO, after eliminating bright short bursts, and compared 
it with that in Paper I. We fixed the spectral parameters at 
the values reported in Paper I, and only three parameters 
of the XISO response (a normalization factor, a slope and 
an offset of the gain) were left free to be fitted. Then, we 
determined the normalization factor, the slope and the gain 
to be 1.079, 0.90, —0.05, respectively. In this paper, we use 
these correction factors. 



APPENDIX B: THE STATISTICAL 
SIGNIFICANCES OF DETECTED BURSTS 

We have already evaluated the detection significances in §2.2 
using only the HXD-PIN data. Here we further evaluated the 
statistical significance of these detections combing the XIS, 
HXD-PIN, and HXD-GSO together. Especially, the XIS in- 
struments are placed at different location in the Suzaku 
satellite, temporal coincidences of the X-ray events between 
the HXD and the XISO make the detection much plausible. 
A possibility to detect n events, from the Poisson distribu- 
tion of its average A count sec^ 1 , during a time duration of 
T can be written as P(A,n) = (XT) n exp(— AT) /nl. Thus, 
multiplied by the independent instruments, the total proba- 
bility of these bursts become Ptotal = Pxis -Ppin -Pgso- Here 
we employ T = 1 sec and corresponding burst rates listed in 
Table [T] together with average rate of Axiso = 3.98 counts 
s _ , Apin = 0.76 counts s~ , Agso = 13.38 counts s" 1 (see 
§2.1). Considering the total bin number of ribm = 4.5 x 10 4 , 
chance occurrence ribmPtotai of these bursts are estimated 
to be below 10 -5 . Therefore, we confirmed these events are 
indeed short bursts. 
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Figure Al. (left) Accumulated light curves of 5 short bursts (Burst-4,7, 9, 14, 16) from XISO (top) and for HXD-PIN (bottom) both 
in the 10-14 keV energy range, (right) Same as the left panels, but from the HXD-PIN (top) and HXD-GSO (bottom) both in the 40-70 
keV energy range. Burst-2, 6, 7, 9, 10, 11, 14, 16 and 18 are used for the accumulation. 



i i i I i i i i I i i i i I i i i i I i i i i I i i i 




Time (ms) 

Figure A2. Cross-correlation coefficients calculated through equation ((AT}. Panel (a) is obtained from the P-sum mode XISO data to 
the HXD-PIN data both in the 10-14 keV energy range, corresponding to the data in figure [ATT left). Panel (b) is obtained from the 
HXD-GSO data to the HXD-PIN data both in the 40-70 keV energy range, corresponding to those in figure [ATT right) . 



APPENDIX C: CHECK OF THE DATA LOSS 

The burst data (except for Burst-6 and Burst- 10), were ex- 
amined for possible losses or pile-ups, through the same 
procedure as applied to the much brighter burst from 
SGR 0501+4516 recorded with Suzaku in 2008 (jEnoto et all 
120091 ). There are three possible forms of the data loss; (1) 
dead time in the event handling at the analog electronics 
(HXD-AE), (2) the limited event transfer rate from HXD- 
AE to the digital electronics (HXD-DE), and (3) that from 
HXD-DE to the spacecraft data processor. Since the peak 



count rates of these bursts are less than ~60 counts (15.6 
ms) -1 , the effect of (1) is estimated to be at most <8%. and 
on average ~1%, since the data acquisition time (20 fxs) is 
much shorter than the burst event interval. Pileup is much 
less effective due to a short time window for the pulse height 
latch. These event data were not lost in the above (2) and 
(3) processes as well, since the event rates were much smaller 
than the maximum limits; 4 kHz for (2) and 400 counts (4 
sec) -1 for (3). Thus, our samples are free from the data loss 
due to the dead-time or the event pile-ups. Since the dead- 
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time effect is typically ~ 1% for most cases, we does not 
correct the dead-time. 



